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ABSTRACT 



A new three-dimensional model of the solar photosphere is presented in this paper and made publicly available to the community. 
This model has the peculiarity that it has been obtained by inverting spectro-polarimetric observations, rather than from numerical 
radiation hydrodynamical simulations. The data used here are from the spectro-polarimeter onboard the Hinode satellite, which 
routinely delivers Stokes /, Q, U and V profiles in the 6302 A spectral region with excellent quality, stability and spatial resolution 
(approximately 0.3"). With such spatial resolution, the major granular components are well resolved, which implies that the derived 
model needs no micro- or macro-turbulence to properly fit the widths of the observed spectral lines. Not only this model fits the 
observed data used for its construction, but it can also fit previous solar atlas observations satisfactorily. 
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1. Introduction 

The controversial debate sparked in recent years regarding the 
chemical abundances of the Sun and other stars, and the pro- 
found consequences that the proposed revision would have in 
many areas of astrophysics, has bluntly put forward an alarm- 
ing weakness in the very foundations of our field. In particular, 
this issue has revealed the pivotal importance of the choice of 
a suitable model atmosphere in the process of abundance deter- 
mination. Traditional one-dimensional (ID) models of the solar 
atmosphere were produced empirically by adjusting their param- 
eters to reproduce more or less adequately a number of observ- 
ables, including spectral lines and continua. Some of the most 
widely used empirical ID models ar e the Harvard-Smithso nian 
Referenc e At mosphere (HSRA, see Ginge rich et al. 197lh. o r 
those of iHoiw eger & MueUei (1 19741) . IVernazza et al.1 (I l98ll) . 
iFontenlaet arj(il993l) . etc. 



A new generation of three-dimensional (3D) models, de- 
veloped theoretically ab initio from radiation hydrodynami- 
cal simulations was used in a series of papers to propose 
a revision of the solar chemical comp osition lea ding to a 
significantly lower metalicity (see, e.g., lGrevesse"et al] l2007l 
and references therein). Of particular importa nce has been 
the issue of the oxygen abundance raised by lAsplund et alj 
(1200 4). prompting what so me authors have dubbed the so- 
lar oxygen crisis dAvres et a l. 2006) and followed by contro- 
versy on whether the proposed revision should be adopted or 
(e.g., 'Bahcall et al.' '2005^; 'Socas-Navarro & Nortoij I " 



not 

Cent eno & Socas-Navarro 2008; Basu & Antia 2008; 
2008t iScott et alj i2009). To add even more confusion, another 



theoretical 3D model was used bv Caffau et al. (2008) to derive 
th e solar oxygen abunda nce, resulting in higher values than those 
oflAsplund et al.' ('2004'). 

^yres et al. (2006) points out that the relative merits by 
which one measures the success of a 3D theoretical model and 
those of an empirical ID one are certainly different. They sug- 



gest that while 3D is obviously to be preferred over ID, an em- 
pirical model is more suited than a theoretical one for abundance 
determinations and it is not clear how these two factors balance 
out in the trade off. 

The path taken in this work is an attempt to combine the best 
of both approaches by deriving a 3D model from observations. 
This is possible now for two reasons mainly: l)we have instru- 
mentation with the capability of providing detailed spectra with 
sufficient spatial resolution to resove the granular motions in the 
solar photosphere; and 2)because the analysis techniques (e.g., 
inversions) are mature enough that it is computationally afford- 
able to undertake a project that involves the detailed study of a 
very large set of profiles. 

Nonwhitstanding the focus of the present discussion on 
chemical abundance determinations, an empirical 3D model is 
also of potential usefulness in a wide range of investigations. 
Noteworthy examples are aiding in the development and veri- 
fication of theoretical ab initio models, predicting the expected 
shapes of continua and lines and, since in this case it also in- 
cludes the photospheric magnetic field, their Zeeman polariza- 
tion signatures as well. 

The model presented here is publicly available and may be 
downloaded both as an IDL savefile or in raw binary formafl 
from the following URL: 

\protect\vrule widthOpt\protect\href {f tp : //download : dal 

The files are licensed under the GPLv3 general public licens^ 
which explicitly grants permission to copy, modify (with proper 
credit to the original source and explanation of the modifica- 
tions) and redistribute the software. 



' As a courtesy, potential users are kindly requested to contact the 
author explaining the nature of their investigations and the intended use 
of the model. 

- See ht tp://www.gnu.org/licenses/gpl.htmll 
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2. Observations and data reduction 

The dataset used in this work was acquired starting at UT 
19:32:10 on 2007 September 24 with the spectro-polarimeter 
(SP) of the solar optical telescope (SOT) onboard the Hinode 
satell i te iKosuei et al.l | 2007|: llchimoto et"all2008l: IShiniizu et aP 



l2008t ISuematsu et al.ll2008l: iTsuneta et al.ll2008 l). The observed 
field of view was very close to disk center, at heliocentric coor- 
dinates (-16, -6) arc-seconds. The slit stepping projected on the 
solar disk was of 0.15" and its length of 162" with a cadence 
of 13 s (exposure time was 12.8 s). A total of 1024 slit positions 
were recorded, resulting in a total field of view of 15 T'x 162" (al- 
though only a smaller sub-field will be analyzed in detail here, 
as explained in Sect.[3]below). 

As usual, the spectral coverage of the Hinode SP has 112 
wavelength samples with a pixel sampling of 2 1 .4 mA and span- 
ning the range between 6300.89 and 6303.27 A. The absolute 
wavelength calibration was obtained by comparing the average 
spectrum to the Kitt Peak Fo urier Transform Spectrometer (FTS) 
disk center intensity atlas of Neckel & Labs ( 1984). With this ex- 
posure time, the signal-to-noise ratio in the data, measured as the 
standard deviation of the Stokes Q, U and V profiles in the con- 
tinuum, is of 1 .2x10"^ in units of the continuum intensity. All of 
the Stokes parameters exhibit approximately the same amount of 
noise. 

All the data were first processed from Level to Level 1 us- 
ing the standard Hinode SOT/SP data reduction pipeline. The re- 
sulting Level 1 data were then subject to a number of additional 
steps, as follows: 

- Remove wavelength offset in the slit stepping direction: The 
absolute wavelength calibration was obtained for the aver- 
age spectrum. However, if one computes the line minimum 
position for each (x, y) point (with x being the slit stepping 
direction and y the direction along the slit), average along 
the y-direction and examine its variation with x, a fluctuating 
pattern appears. The y-averaged line center varies by as much 
as 1 1 pixels (from minimum to maximum) in the x direction. 
This variation is measured consistently in both Fe i lines at 
6301.5 and 6302.5 A (the same value is obtained for both 
lines, with the maximum diff'erence over the whole range be- 
ing 0.02 pixels). To remove this undesirable effect, the line 
center displacement averaged over the two lines and the y- 
direction is subtracted by reinterpolating the four Stokes pro- 
files at each spatial position. As mentioned above this is done 
over the full original field of view of 1 024x 1 024 pixels to en- 
sure that we have as much statistics as possible in the spatial 
average along the slit. 

- Remove pointing jumps in the slit direction: Visual inspec- 
tion of the continuum maps shows some locations with ob- 
vious local discontinuities, caused by the slit moving a few 
pixels up or down in the y-direction. This glitch occurs rela- 
tively infrequently but has been corrected by comparing each 
pixel with its surroundings in the x-direction (3 pixels to the 
left and 3 pixels to the right) and testing if shifting the slit by 
up to 10 pixels up or down would improve continuity in the 
continuum map. The optimal continuity shift found in this 
manner is then averaged along the y-direction (again taking 
the full 1024 pixels to maximize statistics) and applied to the 
data. The results of this procedure may be seen in the exam- 
ple shown in Fig.[T] 

- Following lDanilovic et al.l (l2008l) . a parasitic stray light con- 
tamination of 5% is removed from all the intensity spectra 





Fig. 1. Detail of the continuum intensity in the observed field of 
view showing some glitches due to the sudden displacement of 
the frame along the slit direction (the arrows point to some of 
the most obvious examples). Left panel: Original image. Right 
panel: Corrected as explained in the text. 



before normalizing all the dataset to the average continuum 
intensity over the entire field of view. 

After all data processing, a few flatfield artifacts still remain 
visible in the continuum images. However, the smaller subfield 
selected for further analysis below has been selected to avoid 
the affected areas. The resulting rms granulation contrast in the 
region is 7.7% of the mean continuum value. 



3. Analysis 

For practical reasons, a smaller subfield of view of 200x200 
pixels has been selected for further analysis. This size is large 
enough to include abundant statistics on granulation and quiet 
Sun fields (including both network and internetwork) and at the 
same time small enough to allow a full inversion of each spa- 
tial point individually in a reasonable time span using a super- 
computer Figure |2] shows continuum intensity and a synthetic 
magnetogram (constructed simply as the absolute value of the 
circular polarization signal integrated over a 100 mA range in 
the blue lobe of both lines) of the original full field of view and 
the subfield selected for inversion. 

The inversions were c arried out with the code NICOLE 
(ISocas-Navarro et aT 2010), which is an i mproved implementa- 
tion of the original NLTE inversion code o fH ocas-Navarro et al.l 
dloool) . It makes use of response functions similarly t o the p opu- 
lar LTE code SIR of iRuiz Cobo & del Toro Iniestal H992}, but 
has some features that are important for our purposes here. 
NICOLE is supported by a very wide variety of platforms and 
has built-in MPl parallelization, which makes the use of super- 
computers rather simple and straightforward. Moreover, it has 
NLTE capabilities, including an option that is used for part of 
this work (see Sect. l3.1b in which it is possible to supply a (fixed) 
set of departure coefficients for the upper and lower levels of the 
spectral lines. 

The calculations presented in this paper have been performed 
using the LaPalma supercomputer of the Instituto de Astrofisica 
de Canarias which harbors 512 64-bit processors, running at a 
speed of 2.2 GHz. A typical run, including the various passes 
described below, requires about 6 hours on 50 CPUs. 
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Arc-seconds Arc-seconds 

Fig. 2. Full field of view (upper panels) and subfield selected for inversions (lower panels). Left: Continuum image. Right: Synthetic 
magnetogram, saturated at 5% of the average continuum intensity. The selected subfield avoids some flatfield artifacts that are 
visible as horizontal segments in the upper left panel. 



The observed spectral region contains two prominent Fei 
lines (at 6301.5 and 6302.5 A) whose profile shapes are fitted 
by the inversion code to find the optimal depth stratification of 
the atmospheric variables. These lines are well described with 
a LTE treatment (the effects of NLTE corrections will be ana- 
lyzed in Sect. l3.1l bel ow). The atomic param eters are taken from 
the VAL-D database jPiskunov et alj[l995l) . except for the tran- 
sition probabilities. For the 6301.5 A line we use the laboratory 
measurements of Bard et al. (1991). 

Unfortunately, no similar measurements exist for the 
6302.5 A line and we derived it empirically in the following 
manner We first inverted the 6301 .5 A line in the FTS atlas men- 
tioned above. With the model thus obtained we synthesized the 
6302.5 A line adjusting the log{gf) parameter until a satisfac- 
tory fit to the observations was attained. CoUisiona l broadening 
is treated with the formalism of lAnstee & O'Mara ( 199 3 with 
the par ameters a and cr obtained using the code of B arklem et al.l 
(Il998h . The atomic parameters used are given in Table [T|(ro de- 
notes de Bohr radius). 

For each pixel in the 200x200 subfield, the inversion pro- 
cedure takes a starting guess model atmosphere and iteratively 
modifies it by adding a correction (which, in general, is depth- 



Table 1. Spectral line data 





Excitation 


log(gf) 


Term Term cr 


a 




pot. (eV) 




{lower) {upper) {r^) 




6301.5 


3.654 


-0.718 


^P2 ^£>2 834.4 


0.243 


6302.5 


3.686 


-1.160 


^Pi ^Do 850.2 


0.239 



dependent) to it, seeking the best fit to the observations with the 
synthetic profiles computed from that modified model. The cor- 
rection that is applied to the guess model at each iteration may be 
constant with height (depth-independent), linear with log(T5()(K)) 
(the logarithm of the continuum optical depth at 5000 A) or it can 
be constructed as a concatenation of linear segments spanning 
the whole depth range. In this last case, the number of segments 
to use may be set arbitrarily by the user In this manner, depend- 
ing on the amount of information available, we can decide how 
much detail of the depth stratification we wish to retrieve. One 
needs to reach a compromise to adjust the amount of freedom 
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in the model to the information available in the data. Too much 
freedom results in degeneracies and leads to uniqueness issues 
and other complications. On the other hand, being too restrictive 
results in poor fits and not extracting all of the available infor- 
mation. Usually, some experimentation with a few test cases is 
very useful to determine a nearly optimal set of parameters. 

The choice of freedom in each physical parameter of 
the model (temperature, magnetic field, velocities, etc) is 
made in practice by selecting the ri umber of inversion nodes 
dRuiz Cobo & del Toro Iniesta Il992h . which are actually the 
free parameters in the code. Using one node we have a 
depth-independent correction for the entire atmosphere. With 
two nodes we produce a correction that scales linearly with 
log(T5ooo). Similarly, with three or more nodes we produce 
more complicated variations . Following the recommendation of 
[Ruiz Cobo & del T oro Inies ta ( 1992), we proceed in two succes- 
sive cycles to improve convergence. We start in the first cycle 
with a relatively small number of nodes to obtain a first approx- 
imation to the solution, fitting the overall shape of the profiles. 
We then increase the number of nodes to obtain a better solution 
allowing more freedom to fit more subtle properties, such as line 
asymmetries. 

In addition to these cycles, NICOLE implements a multiple- 
initialization scheme in which the inversion (both cyles) is re- 
peated a number of times randomizing the initial guess. The 
process is stopped when a good fit is achieved or when a preset 
maximum number of inversions have been done. The criterion to 
decide what we consider a good fit and the maximum number of 
inversions are both defined by the user This is helpful when, as 
in this case, one is batch processing a large number of profiles to 
minimize the risk of the algorithm settling in a secondary mini- 
mum. The price to pay for the ensuing improvement in stability 
is, of course, more CPU time. In the calculations presented here, 
the stopping criteria are: a)a fit better than 1% of the profile on 
average or b)up to five inversion attempts. 

As a starting guess we take the HSRA model. Two different 
approaches are taken depending on whether the profile to invert 
exhibits significant polarization or not. In the first case, which 
accounts for 34% of the profiles, we consider two atmospheres 
coexisting side-by-side in the pixel, one magnetic and the other 
non-magnetic. To determine the non-magnetic component, we 
start with a preliminary inversion of Stokes / only (weights for 
Q, U and V are set to zero) without any magnetic field. In this 
inversion we allow for some micro-turbulent line broadening to 
improve the fit since it is unlikely that the external atmosphere 
will be entirely field-free. However, we do not invert for a mag- 
netic field here since our main focus is the field in the second 
component which is the one that gives rise to the Stokes Q, U 
and V profiles. The resulting model and Stokes / profile from 
this preliminary inversion are then taken as the (fixed) external 
non-magnetic atmosphere in a subsequent inversion in which the 
filling factor is a free parameter In the case that the pixel con- 
sidered does not exhibit Stokes signal we proceed with a sim- 
ple one-component inversion. The number of nodes employed 
in both cases is shown in Table |2] As mentioned above, only the 
first (preliminary) inversion has microturbulence so that the code 
is able to handle the extra broadening due to the magnetic field. 
In the second and third inversions, where the magnetic field is 
taken into account, microturbulence is set to zero. 

In the table, Bio„g refers to the longitudinal (i.e., along the 
line-of-sight) component of the magnetic field, whereas and 
By are the transverse components, projected on the plane of the 
sky. The x-direction is the field azimuth reference, defined by 
the plane of vibration of light with Q>Q and U -Q in the cali- 



bration pipeline. The synthetic Stokes spectra computed at each 
step of the iteration are convolved with the instrumental profile 
of the Hinode SP (Lites, private communication) to make sure 
that synthetic and observed profiles are comparable. 

A final pass is performed restarting the inversion with a 
smoothed version of the result to fix a small percentage of pix- 
els that did not converge properly. The smoothing is done with a 
3x3 median box, but excluding the central point (i.e., the pixel 
to be re-inverted). 

A sample of the resulting cubes can be seen in Fig.|3]which 
shows horizontal (in tsooo) cuts of the temperature and line-of- 
sight velocity at three representative heights in the photosphere. 
For the two-component inversions (pixels with magnetic signal), 
an average of the internal and external atmospheres weighted 
with their respective filling factors is shown. 

The figure shows the characteristic convective motions at the 
base of the photosphere (top panels), with hot granules upflow- 
ing and cool lanes downflowing. The flow structure does not 
change too much until we reach the upper photosphere (bottom 
panel) in which the granulation pattern starts to dissolve leaving 
structures that vaguely appear to be vertically oriented. Since 
there is no preferred direction in the field of view, such struc- 
tures are likely the result of observing a very dynamic upper 
atmosphere through a vertical slit. The temperature stratification 
starts with the granulation at the bottom but at log(T5ooo)= -1 
is dominated by hot patches in the magnetic areas embedded 
in a more uniform non-magnetic background. Finally, at the 
log(T5ooo)= -2 height we start to see reversed granulation, where 
the photospheric granules are now cooler and the lanes are hot- 
ter. 



3. 1 . NLTE corrections 

The Fe i lines used in this work are affected by (generally very 
small) NLTE efi'ects (Shchukina & Trujillo Bueno 2001). Both 
the upper and lower levels of the transitions are slightly under- 
populated with respect to their LTE values. This has virtually no 
effect on the source function, which goes with the ratio of the 
level populations, but produces a small overall decrease in the 
line opacity. 

We introduce a correction for NLTE eflFects in the inver- 
sion sch eme by using some departu re coefficients computed 
by Shchukina & Trujillo B uenol (1200 II) in the 3D hyd rodynamic 
model of Asplund et al. (2003 see also lStein & Nord lund 1998i). 
We take the run with log(r50()o) of the departure coefficients for 
both levels of the 6301.5 and 6302. 5 A transitions computed by 
IShchukina & Truiillo Buenol (1200 ih in a typical granule and a 
typical lane in the 3D simulation. These values are then assigned 
to a particular granule and a lane pixel in the Hinode obser- 
vations. For all other pixels we carry out a linear interpolation 
based on the continuum intensity. The departure coefficients ob- 
tained are kept fixed during the inversion, even though the at- 
mospheric parameters change. While this method is not exact, it 
is a good approximation to correct for an effect that is already 
small anyway. A full NLTE solution for the Fe atom would be a 
tremendous undertaking and a full 3D inversion such as the one 
presented here would not be possible in practice due to the pro- 
cessing power required. With our approximation, the departure 
coefficients may be inconsistent with the starting guess atmo- 
sphere but as the inversion proceeds and approaches the solution, 
they become gradually more realistic. 




Fig. 3. Horizontal cuts of the temperature (left) and line-of-sight velocity (right) at log(T5ooo) - (top), -1 (middle) and -2 
(bottom). Units are kK for temperature and km s"' for velocity. The astrophysical convention, where positive velocities indicate 
redshift (downflows), is adopted throughout this paper. 
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Table 2. Inversion nodes 







Magnetic pixel 






Non - magnetic pixel 




Physical 


Inversion 1 




Inversion! 








Parameter 


Cycle 1 


Cycle! 


Cycle I 


Cycle! 


Cycle 1 


Cycle! 


T emperature 


3 


5 


2 


3 


3 


5 


L.o.s. velocity 


1 


3 


1 


2 


2 


3 


Microturbidence 


1 


1 














^tong 








1 


2 
















1 


2 








By 








1 


2 








Filling factor 








1 


1 









4. Results 

With the NLTE correction, the cores of the Hnes are fitted better 
and therefore the results for the upper layers are less noisy (see 
Fig. m compare bottom-left panel to that in Fig. [S). Figure |4] 
shows also the magnetic flux density in the region, revealing 
many internetwork field structures. The transverse flux density 
is more noisy because the Stokes Q and U signals are weaker 
than Stokes V and correspondingly more difficult to detect. 
Nevertheless, we can see that a large fraction of the area harbors 
fields with a measureable horizon tal component, in agreement 
with recent results such as those o f iLites et al.l ( l2008h . 

The magnetic field height gradient is mostly negative as one 
would expect. However, some concentratios show an increase 
of the field strength with height. Temperatures inside the mag- 
netic atmospheres always decrease with height but there is a 
wide range of variation of ~1 kK among various pixels. The 
mass flows may be upward or downward directed. However, 
the stronger flux elements are associated with downflows. If we 
consider pixels with a longitudinal magnetic flux stronger than 
100 G, 84% harbor downflows and only 16% are upflowing. The 
external (non-magnetic) atmosphere is also undergoing down- 
flows in those pixels. 

Generally speaking, the fits to the individual profiles are 
good. Figure |5] shows the average intensity profile in the ob- 
served region and the average synthetic profile from the model. 
Notice that this is not actually a fit but an average of the 200x200 
individual fits. 

We can also test the model against atlas observations. If we 
synthesize the profiles in the absence of macroturbulence and 
instrumental profile (i.e., with very high spectral resolution), we 
should obtain a profile very similar to that of an intensity at- 
las. In F ig. |6] we can see the comparison with the FTS atlas 
dNeckel & Labs 1984). Overall, the model reproduces the atlas 
observations with very high fidelity, including the line broad- 
ening. The absence of microturbulence does not represent any 
problem for obtaining realistic line shapes. This is not entirely 
surprising, since the observations have sufficient spatial resolu- 
tion to resolve the granular components (granulation is the main 
contributor to the line broadening traditionally characterized by 
microturbulence in ID models). It is, however, a noteworthy fea- 
ture since the success in reproducing line widths without requir- 
ing any microturbulence has been remarked as a strong argument 
supporting theoretical 3D models derived in recent years. 

When we compare the average temperature stratification of 
our model to others previously existing in the literature, we find 
that it is similar, although slightly cooler in the middle photo- 
sphere (some 100-200 K at log(r5oo())= -1) and warmer at the 
top (~150 K at log(T5()oo)= -3). Figure [T] compares the horizon- 




6302.0 
Wavelength (A) 



Fig. 5. Average intensity spectrum from the Hinode SP observed 
dataset (solid) and from the 3D model (dashed). Spectra are nor- 
malized to the continuum intensity (Ic). 



tally averaged temperature to t hat of the Harvard-Sm ithsonian 
Reference Atmosphere ( HSRA, Gingerich et al. M97Tb and to the 
average temperature of lAsplund et al. (2004). The HSRA is a 
semiempirical ID model and one expects some differences with 
the one presented here because in principle the ID model that 
fits a given spectral dataset is not necessarily the same as the av- 
erage stratification of a 3D model that fits the same dataset. In 



other words, the HSRA already inco rporates an ernpirical cor- 
rection for 3D effects. The model of lAsplund et al.l (l2004l) is a 
3D model but obtained from theoretical simulations. Compared 
to HSRA, it also produces an elbow around log(T5(K)o)= -1 but 
not as pronounced as our model. In the upper layers, on the other 
hand, both HSRA and the Asplund et al models are very simi- 
lar in spite of the ID vs 3D difference in approach and both are 
cooler than the present model. 

The importa nce of 3D effects ha s been emphasized in re- 
cent work (e.g., lAsplund et al.l 120041) in the context of chemi- 
cal abundance determinations. To assess the importance of 3D 
effects in our model we have compared the average emerging 
spectrum from what would be obtained if one first averages the 
atmospheric stratification and then compute the spectrum of such 
ID average model. The results, displayed in Figure |8] are in this 
case rather small (but noticeable), with the main departure been 
at the core of the lines (i.e., affecting the upper atmospheric lay- 
ers). The lines computed in ID are also slightly narrower, which 
goes in the direction of compensating the change in equivalent 
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Fig. 4. Horizontal (in the T5000 depth scale) cuts of several interesting physical parameters in the multi-cube resulting from the 
inversion, including NLTE correction. The left column shows the temperature in kK at three different heights. The top row also shows 
the line-of-sight magnetic flux density (middle panel) in G and the transverse (on the plane of the sky) magnetic flux density (right 
panel) at the base of the photosphere. The middle row shows the variation of the flux density from log(Tg(K)o)= to log(T5ooo)= -1 
(middle panel: line-of-sight component; right panel: transverse component). The bottom row shows the variation of the temperature 
in kK (middle panel) and the line-of-sight velocity in km s ' (right panel) from log(T5ooo)= to log(T5ooo)= -1 inside the magnetic 
atmosphere. 



width. In fact, the difl'erence in equivalent width between the ID 
and 3D profiles is at the 1.5% level. 



5. Uncertainties 

It is important to estimate the errors and uncertainties in the 
derived model. Systematic errors may come from the normal- 
ization of the Hinode data since the observations are not pho- 
tometrically calibrated. Instead, relative intensities to the disk- 
center average continuum have been employed, with the sub- 
sequent assumption that such disk-center average intensity is 



well reproduced by the HSRA model. The visible continuum 
is very sensitive to temperature variations, which has the ad- 
vantage that it is relatively easy to determine the photospheric 
temperature. On the downside, a small error in the intensity 
calibration might have relatively large eff'ects on the tempera- 
ture determination. For instance, a 1% error in the continuum 
calibration results in a AT of between 150 and 200 K (it is 
slightly different for granules and lanes). This example is par- 
ticulary pessimistic, as we do not expect such a large error in 
the average intensity determination. Furthermore, the average 
temperature of our model at tsooo = 1 (6441 K) is consistent 
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Fig. 6. FTS atlas intensity spectrum (solid) and average profile 
synthesized from the 3D model (dashed) without instrumental 
profile (i.e., with very high spectral resolution). The lines at 
6302.00 and 6302.76 A are telluric lines (absent in the Hinode 
data). Spectra are normalized to the continuum intensity (7^). 




Og(T500ll) 



Fig. 7. Comparison of the average temperature stratification in 
our 3D model (solid l ine) to HSRA (dash ed) and the aver- 
age stratification of the lAsplund et alj (|2004|) 3D model (dash- 
dotted). 



with values from several other models existing in the litera- 
ture, such as HSRA (6390 K) o r thos e of [Asplund et al.' (2004', 
6412 K ), Holweger & Muefle^ (flOTl 65 30 K), .Vernazza et ah 



(I1981L 6424 K1 and lpontenla et al.l (I1993L 6520 K). 



5. 1 . Sensitivity range 

Semi-empirical models derived from inversion of spectro- 
scopic data, such as this one, are reliable only in the atmo- 
spheric region where the spectral lines exhibit some sensi- 
tivity to the physical parameters. A suitable way to explore 
this range is by using the so-calle d response functions (e.g., 
iRuiz Cobo & del Toro Iniestalll994l) . Mathematically speaking, 
the response function of a given Stokes parameter at a given 



Fig. 8. Comparison of the average spectrum from the 3D model 
(solid) to the spectrum produced by the average model (dashed). 




0.0 0.5 
A - 6302 (A) 




0.0 0.5 

A - 6302 (A) 



Fig. 9. Normalized response functions R(A, r) to temperature in 
a granule (upper panel) and a lane (lower panel). 



wavelength 1(A) to perturbations in a given physical parameter 
at a certain height in the atmosphere T(t) is expressed as: 



R{A,T) 



dI(A) 
dT{T) 



(1) 



Here we use the same nomenclature to refer to the discrete form 
(replacing the derivative with a ratio of finite increments) of 
Eq ([T]) which, albeit not entirely equivalent from the mathe- 
matical point of view, is very convenient in practical terms for 
our study. The most straightforward way to compute a response 
function is by brute force from its definition. We start with a ref- 
erence atmosphere from which the emerging spectrum is known 
and then perturb this atmosphere, one point at a time, recompute 
the spectrum and calculate the ratio M(A)I Is.T{t) (assuming that 
the perturbed variable is temperature at depth r). Response func- 
tions to temperature in a granule and an intergranular lane are 
displayed in Fig.|9l 

From the response functions we can already see that the sen- 
sitivity range of these two lines spans from roughly log(r5ooo)= 
to -3, approximately. Let us go a step further and use this infor- 
mation to estimate the actual errors in the inversions. To this 
effect it is important to consider not only how sensitive the lines 
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-3.0 



-2.5 -2.0 -1.5 



Fig. 10. Estimated uncertainty in the temperature stratification 
for three different model atmospheres. SoHd line: Brightest con- 
tinuum location in the field of view (a granule). Dashed line: 
Darkest continuum location in the field of view (a lane). Dotted 
line: A location having a brightness equal to the average of the 
whole region. Coordinates on the lower panel maps of Fig. |2]are 
(14.5,8.2) for the granule, (12.0,25.5) for the lane and (22.6,8.5) 
for the average brightness point. 



are but also how good a fit we have obtained. If we had only one 
wavelength, it is straightforward to see that the error in the de- 
termination of a given parameter (e.g., T to fix ideas) may be ap- 
proximated by Ar - AI/R, where AI is the difference between 
the observed and the synthetic observable and R is the response 
function. 

Consider now a case with many wavelengths and assume 
for simplicity that the inversion result may be regarded as 
a weighted average of individual measurements of ^(t), one 
at each wavelength, all weighted by their individual response 
R(A, t). We can then use the expression for the variance of a 
weighted mean to obtain in our case: 



1 



Ar(T)2 



AI(A) 



(2) 



FigurefTOlshows the uncertainties in the temperature determina- 
tion obtained using this expression in three different pixels of the 
field of view, which have been selected according to their con- 
tinuum brightness to include a granule (bright point), an inter- 
granular lane (dark point) and an average location. The residual 
discrepancy between the average observed and synthetic spectra 
is introduced in the equation as the AI{A) function. The figure 
shows that extremely small errors (~10 K or smaller) are ob- 
tained for a fairly wide range of heights, between log(T5ooo)= 0.5 
and -1.7, approximately. Up until log(T5oo())= -2.5 we can still 
have relatively small errors (~150 K) but above that height, mea- 
surements should be regarded as highly uncertain. Notice, how- 
ever, that the dependence of the errors with the model atmo- 
sphere chosen is significant for the upper layers. The granule 
model is the one that has the error increasing more rapidly with 
height between log(Tgooo)= -2 and -3. In the lane model the 
error increases more slowly. Finally, in the average atmosphere 
one can reach log(r5oo())= -3.4 with very small errors of ~50 K. 



Fig. 11. Models obtained when inverting the synthetic average 
spectrum with different values of the Fe abundance. The solid 
line is for the reference value of 7.50. The dashed lines are for 
7.40 (lower curve) and 7.60 (upper curve). 



5.2. Abundances 

For the calculations in this work we have employed the sola r 
chemical composition published by 'Grev esse & Sauvall (119981) . 
Since the spectral lines analyzed are from Fe transitions, this 
is the abundance value tha t will have a stronger imp act on our 
results. The 7.50 value in iGrevesse & Sauvall (11998 ) seems to 
be very well established and, when NLTE effects are considered 
dShchukina & Truiillo Buenol2001h . in good agreement with the 
meteoritic value. The rest of the elements are relevant only for 
the calculation of the background opacities. 

To assess the impact of possible uncertainties in the abun- 
dances employed, the synthetic profile produced from the aver- 
age model was inverted with different sets of abundances. A first 
test had the Fe abundance fixed and a constant scaling was ap- 
plied to all other elements. The scaling factor ranged from 0.75 
to 1.25. For the second test the Fe abundance was varied be- 
tween 7.40 and 7.60, while the rest of the chemical composition 
was kept constant. 

As expected, the only significant difference was found in the 
test where the Fe abundance was modified. Figure [TT] shows the 
inversion with the reference value of 7.50 compared to the mod- 
els obtained with the interval extremes (7.40 and 7.60, respec- 
tively). The changes are very small, mostly because the emerg- 
ing intensity is extremely sensitive to temperature in this range. 
This means that a very small variation in the model tempera- 
ture is able to compensate for significant variations in the line 
parameters. The upside of this conclusion is that semiempiri- 
cal models such as the one presented here are rather robust and 
not very sensitive to uncertainties in line data, abundances, etc. 
Unfortunately, the downside is that abundance determinations 
are extremely sensitive to the model employed. 



6. Conclusions 

The model presented here is novel in that it is three-dimensional 
in nature and has been obtained from real observations. It may 
hopefully serve as another platform for studies of chemical com- 
position (in which more 3D models are urgently needed), to pro- 
duce synthetic spectra (with or without polarization), analyze 
properties of the quiet Sun and their height variation or to learn 
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about the physics of (magneto-)convection from direct compari- 
son with numerical simulations. 

This work may be improved upon by using new observa- 
tions with even higher spatial resolution or more spectral fea- 
tures. Existing inversion codes such as NICOLE have the capa- 
bility of incorporating lines from different elements and even to 
combine LTE and NLTE lines. It should be possible, in principle, 
to combine photospheric and chromospheric lines with the aim 
of producing a 3D model like this but spanning a much greater 
range of heights, all the way up to the middle chromosphere. The 
main limitation, however, is that we currently lack adequate in- 
strumentation to produce such high-resolution multi-wavelength 
observatio ns. Future planned ground-bas ed telescopes such as 
the ATST (iKeller et al.ll2002h or the EST dCollados et al.ll20Toh 
or space-born observatories such as the plan B option of Solar-C 
(Hinode's successor, still in early planning stages) would pro- 
vide an enormous leap in our ability to acquire the necessary 
data. 

Acknowledgements. The author thanks Nataliya Shchukina and Javier Trujillo 
Bueno for kindly providing the departure coefficients used for the NLTE corTec- 
tion in this work. 

Hinode is a Japanese mission developed and launched by ISAS/JAXA, col- 
laborating with NAOJ as a domestic partner, NASA and STFC (UK) as inter- 
national partners. Scientific operation of the Hinode mission is conducted by 
the Hinode science team organized at ISAS/JAXA. This team mainly consists 
of scientists from institutes in the partner countries. Support for the post-launch 
operation is provided by JAXA and NAOJ (Japan), STFC (U.K.), NASA, ESA, 
and NSC (Norway). 

The author thankfully acknowledges the technical expertise and assis- 
tance provided by the Spanish Supercomputing Network (Red Espaola de 
Supercomputacin), as well as the computer resources used: the LaPalma 
Supercomputer, located at the Instituto de Astrofsica de Canarias. 

Financial support by the Spanish Ministry of Science and Innovation 
through project AYA2010-18029 (Solar Magnetism and Astrophysical 
Spectropolarimetry) is gratefully acknowledged. 



Shimizu, T, Nagata, S., Tsuneta, S., et al. 2008, Sol. Phys., 249, 221 
Socas-Navarro, H., de la Cruz, J., Asensio Ramos, A., Trujillo Bueno, J., & Ruiz 

Cobo, B. 2010, A&A, in preparation 
Socas-Navarro, H. & Norton, A. A. 2007, ApJ, 660, L153 
Socas-Navarro, H., Trujillo Bueno, J., & Ruiz Cobo, B. 2000, ApJ, 530, 977 
Stein, R. F & Nordlund, A. 1998, ApJ, 499, 914 

Suematsu, Y., Tsuneta, S., Ichimoto, K., et al. 2008, Sol. Phys., 249, 197 
Tsuneta, S., Ichimoto, K., Katsukawa, ¥., et al. 2008, Sol. Phys., 249, 167 
Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635 



References 

Anstee, S. D. & O'Mara, B. J. 1995, MNRAS, 276, 859 

Asplund, M., Grevesse, N., Sauval, A. J., AUende Prieto, C, & Kiselman, D. 

2004, A&A, 417, 751 
Asplund, M., Ludwig, H., Nordlund, A., & Stein, R. F 2000, A&A, 359, 669 
Ayres, T. R. 2008, ApJ, 686, 731 

Ayres, T. R., Plymate, C, & Keller, C. U. 2006, ApJS, 165, 618 
Bahcaff, J. N., Basu, S., Pinsonneault, M., & SereneUi, A. M. 2005, ApJ, 618, 
1049 

Bard, A., Kock, A., & Kock, M. 1991, A&A, 248, 315 

Barklem, P S., Anstee, S. D., & O'Mara, B. J. 1998, PASA, 15, 336 

Basu, S. & Antia, H. M. 2008, Phys. Rep., 457, 217 

Caffau, E., Ludwig, H., Steffen, M., et al. 2008, A&A, 488, 1031 

Centeno, R. & Socas-Navarro, H. 2008, ApJ, 682, L61 

Collados, M., Bettonvil, F, Cavalier, L., et al. 2010, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. in press. Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference Series 
Danilovic, S., Gandorfer, A., Lagg, A., et al. 2008, A&A, 484, L17 
Fontenla, J. M., Avrett, E. H., & Loeser, R. 1993, ApJ, 406, 319 
Gingerich, O., Noyes, R. W., Kalkofen, W., & Cuny, Y. 1971, Sol. Phys., 18, 347 
Grevesse, N., Asplund, M., & Sauval, A. J. 2007, Space Sci. Rev., 130, 105 
Grevesse, N. & Sauval, A. J. 1998, Space Science Reviews, 85, 161 
Holweger, H. & Mueller, E. A. 1974, Sol. Phys., 39, 19 
Ichimoto, K., Lites, B., Elmore, D., et al. 2008, Sol. Phys., 249, 233 
Keller, C. U., Rimmele, T. R., Hill, F, et al. 2002, Astronomische Nachrichten, 
323, 294 

Kosugi, T., Matsuzaki, K., Sakao, T., et al. 2007, Sol. Phys., 243, 3 
Lites, B. W., Kubo, M., Socas-Navarro, H., et al. 2008, ApJ, 672, 1237 
Neckel, H. & Labs, D. 1984, Sol. Phys., 90, 205 

Piskunov, N. E., Kupka, F, Ryabchikova, T. A., Weiss, W. W., & Jeff'ery, C. S. 

1995, A&AS, 112, 525 
Ruiz Cobo, B. & del Toro Iniesta, J. C. 1992, ApJ, 398, 375 
Ruiz Cobo, B. & del Toro Iniesta, J. C. 1994, A&A, 283, 129 
Scott, R, Asplund, M., Grevesse, N., & Sauval, A. J. 2009, ApJ, 691, LI 19 
Shchukina, N. & Trujillo Bueno, J. 2001, ApJ, 550, 970 



